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Abstract 

Investigating leptonic decays of 'k~,K~, vr®, iC£, we present new bounds on 
some products of two i?-parity violating coupling constants. For mesons of a similar 
structure but so far poor experimental data we give the corresponding formulae, to 
be used in the future. 


1 Introduction 

The MSSM+flp is obtained from the MSSM by adding the following terms to the super¬ 
potential (c./. ref.[T]) 

= i A.,, LI Li S’> A'y, L' Ql of 

+ i c"*" Kjt Uf Df of + e"” Ki 4 Ht (1) 

H, Q, L represent the leftchiral S'?7(2)Tv-doublet superhelds of the Higgses, the quarks 
and leptons; U, D, E represent the rightchiral superhelds of the u-type quarks, d-type 
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quarks and electron-type leptons, respectively; a superscript denotes charge conjuga¬ 
tion; a,b and x,y,z are SU{2)y/- and ^[/(Sjc-indices, respectively; i,i,k and later also 
f,g,l,n are generational indices (summation over repeated indices is implied); 5^^ is the 
Kronecker symbol, £ ■ symbolizes any tensor that is totally antisymmetric with respect 
to the exchange of any two indices, with = 1. The coupling constants \jklKjk 
antisymmetric with respect to the exchange of the hrst two/last two indices. The last 
term in eq. m can be rotated away utilizing a unitary held-redehnition. 

Good agreement between SM theory and experiment gives stringent upper bounds 
on the extra 45 coupling constants Xijk, X'ijk and A"as well as on products thereof. For 
a list of references and the processes dealt with, see e.g. ref. In particular, in 

there are new operators for leptonic meson decays. The SM theoretical predictions for 
the decay widths of mesons and the measured values match up within the experimental 
uncertainty. We can thus determine yet further tight constraints on several products 
of coupling constants: X*X and A*A. This was hrst done in ref. [5] for single coupling 
constants and later in ref. for some products, however only treating charged pions 
decaying via either d-type squark or slepton exchange, respectively. Ref. [Zj treated general 
leptoquark reactions of several particles, one of them the this result was quoted in 
terms of ^p by ref. jH] ; the same result was reached by ref. jD] . Ref. PI among other things 
dealt with the decay of K^, however with only -u-squark exchange contributing to SM- 
allowed processes. The decays of neutral and charged R-mesons were treated in ref. m 
and ref. [12], respectively. We generalize these calculations, focusing on products of two 
coupling constants, and stress where we obtain new or stricter bounds. 

2 ftp-Decay of Charged Mesons 

2.1 Calculation of the Decay Rate 

Consider a negatively charged meson tt®-’ at rest made of a d-type quark d* and a u-type 
antiquark which decays into an antineutrino z/”*" and a charged lepton t'-f, i.e. 

; qp3)), (2) 

the Pi, 2,3 being four-momenta. We now calculate the partial decay rate of this process. 
Focusing on the Yukawa-couplings of the hrst two terms in eq. m leads to, again with 
summation over repeated indices implied, 

D X,jk (^ Pl p i?* + ¥ PlP PlP 4^ ) 
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Figure 1; The tree-level MSSM+^p processes contributing to the decay of the charged 
mesons. 


+ Ahfc ( Pl d/ + # Pl d^ vP + d^ Pl dp ) 

- Kjk Pit Ip* + dF PlPJP + dP PLtlpi) + c.c. (3) 


All spinors are Dirac spinors, the overbar denotes the Dirac adjoint, Pl,r are the pro¬ 
jection operators on the left-/right-handed parts. The fermions are mass-eigenstates. A 
tilde denotes a scalar; the scalars’ subscripts L, R indicate the chirality of the correspond¬ 
ing Weyl spinor. The term in eq. (H together with the c.c. of the 7*^ term, and the 3^'’* 
term together with the c.c. of the 8*^ term lead to the meson decay processes depicted in 
fig- in which give the effective Hamiltonians 

^ V J! fpY PLd\ 

2 ^ mi — 

^ dp 

nP = - Y' Jf ^ p^ ^ 4 ^ 

, m~ 

k 

where m is the mass of a particle. To obtain the hrst equation we employed a Fierz- 
identity. These two Hamiltonians have to be added to the effective Hamiltonian for the 
SM-process 

^ p Pl ^ Y Pl d\ (5) 

Here Gi? is the Fermi constant and Vji is an element of the CKM-matrix. We obtain for 
the transition amplitude Aiij/n 


Mijfn S'^ipi -P2-P3) 


( ef ; I (n^ + + np) 


1 

2TTi 
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d‘^x. (6) 








We expand the fields in the initial and final state, perform the integrations and use 
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is the meson decay constant^ Thus 
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( 8 ) 


Uf ,V^ are the Fourier coefficient functions of ^ respectively. Next we take the 

absolute value square, average over the spins and use the trace theorems. Then we sum 
over n, because the experiments that measured the partial decay widths did not determine 
the flavour of the antineutrinos,^ resulting in 


m'jf {^fn + 26fnRe[Kijfn] + \Kijfn\'^), ( 9 ) 


with 
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( 10 ) 


dp " ' ' iL- 


containing all ftp contributions; 2 Re[iFjj 7 „] in eq.® is due to the interference between 
SM and ^ amplitudes. For simplicity we neglect the phase of the CKM-matrix. The 
partial decay rate is then 


-p SM-\-^p 


X (l + 2Re|Ay„l + E 
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with being an arbitrary antineutrino, and 


-pSM 

^ TT^j^efuf^ 
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( 11 ) 


( 12 ) 


the correction factor Gijf of 0(1) is due to higher order electroweak leading logarithms, 
short distance QGD corrections, and structure dependent effects, see ref.|T2] and also 
ref. jI3] . 

^There are several ways of defining the meson decay constant, differing by factors of v^; in the 
convention we use /tt = (92.4 ± 0.3) MeV, see ref. ESI 

^The upper experimental bounds on tt —> yPe'^' and K yPe^, see ref. m, come from a different 

type of experiment, compared to the one used to determine the branching ratios for tt~ —> yp'^ and 
K~ yp^ . They do not lead to better bounds on the coupling constants. 
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2.2 Calculation of the Bounds 


We prefer not to compare the experimental data directly with eq. (HU, since f^^ij has quite 
a large error. This leads to very weak bounds on Kijfn- To avoid this, we introduce 




l/C 

T^JJ— I/C 


(13) 


with > rriif. If the experimental and SM-theoretical decay rates agree well we have, 
see eq. (HD), \2Re[Kijff] + En <1- Plugging eq.d!!]) into eq.dHl) one gets 

,SM+ftp 

T. (14) 


K 






1 + e^^ij ~ 1 + 2 Ile[Kijff Kijgg\ + ^ \Kijfn\ 


Let A... symbolize the theoretical or experimental uncertainty. If the theoretical predic¬ 
tion lies within the experimental range TZ^f ±ATZ^f, one has 


min 






- 1 A A 


qp exp / ^ exp ' 
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max 


(15) 


We could use this to determine a bound on this general combination of |?p coupling 
constants; however, the bounds on individual coupling constants are typically of the 
order (P( 10 “^), see ref.j 2 ], and thus we limit ourselves to at most two non-zero coupling 
constants at a time, and in each case suppose the other 34 A, A' coupling constants vanish 
feq. (pT)|l . eq. lfTTj) and eq. lfTHjl are also valid for / ^ with -i-i- — 


e™" < 2 Re[iL,,;;] + and forn 7 ^ / < e^. (16) 


We assume that the imaginary parts of the coupling constants are approximately the 
same as the corresponding real parts.^ With Gp = (0.116639±0.000001) x (100 GeV)“^, 
see ref. HSI, we obtain 

-0.330 lEil (^1 ± 2 e™f ± 1 ) < 

^'fik] -2 m±, Re[X',*g, Xfkf] 

(m~/100 GeV)2 ’ m^f {myj + nidi) (m^~/100 GeV)^ 

< 0.330 lEil (^1 ± 2 e-r - 1 ) (17) 


and for n 7 ^ / 

l^'fjk ^nik\ 


2 m 




l^kji ^nkf\ 


{m —j/100 GeV)^ ’ rriif {mui + nidi) (m~^/100 GeV)^ 


< 0.66 \Vj 


3i\ 


'7r*J 


(18) 


' ..... 

The prefactor 2nil.ij/[ni£f{niy^j ± nidi)] results in much tighter bounds for A)))j Xnkf- We 
will apply these results only to processes with sufficiently small experimental error bars. 
^If the imaginary part vanishes the bounds are weaker by a factor of 0(1). 
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2.3 TT 


£/,5 + 
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As a first application, we consider pion decay with /, i, j = 1, g = 2. The SM gives the 2cr 
theoretical value = (1.2354 ±0.0004) x 10“^ (see ref.jHj; the uncertainty mainly de¬ 
rives from Cm and C 112 ). From the partial decay widths at the 2a level in ref. na, namely 

(1.230±0.008) x IO-* and r“iL = 0.9998770 ±8 x 10 


-7 


one calculates = (1.230±0.008) x lO'h Hence, = -0.0107 and = 0.0022. 
With iHiil = 0.9750 ± 0.0008, ref.jS] obtained bounds on a single coupling constant; 
this was updated in ref. uni- We have reproduced their results. The experimental data 
have only marginally changed and the new bounds are < 0.027 m —j/100 GeV 

dR 

and |A 2 ifc| < 0.059 m —r/100 GeV. We obtain bounds for the products of couplings 

dR 




< 0.03 (m~/100 GeV)2, A(jiJ < 0.03 (m~/100 GeV)^ and A(jiJ < 
0.066 (m—^/lOO GeV)^. The hrst bound is redundant since the product of the single 

dR 

bounds is stronger; the second and the third bound are almost the same as the sin¬ 
gle bound on and |A 2 ifc|. Furthermore, we obtain the following new bounds us¬ 

ing me = (0.510998902 ± 2.1 x 10“®) MeV, m^ = (105.6583568 ± 5.2 x lO”®) MeV, 
m^- = (139.57018 ± 0.00035) MeV and mu + md = (8.5 ± 3.5) MeV,^ see ref.[II 
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(19) 


The upper bound/bounds we obtained for Re[A'3\^ A232]/Re[A'{\]^ A212], A312I are 

weaker than the products of the two bounds on the single coupling constants, see ref. |1]; 
there also much stricter bounds were stated for |A).\^ Aifc2| as well as for lA'g^j^ A231I. 


2.4 K- + jyC 

Next we consider charged kaon decay with f,j = 1, g,i = 2. According to ref. jl4j. 
= (2.472 ± 0.002) x 10“® at the 2a level. Experimentally /^k- total ~ 

"^This and rrig are the biggest sources of inaccuracy, going linearly into the bounds on A'* A. 
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and /^^K- total ~ 0.6351 ± 0.0036, at the 2a level jT 


X 10 ° and e™” = —0.10 and = 0.076. Using 


(1.55 ± 0.14) X 10-® 

Therefore = (2.44 ± 0.22) 

IU 12 I = 0.222 ± 0.004 we obtain \\'*^^ A'gafel < 0.04 (m~/100 GeV)2, < 

dR 

0.046 (m^/100 GeV)^. As for the pion these bonnds are almost the same as the ones 
on lAgifc]. Onr bonnds on Re[A'(\^ '^ 12 ^]’ -^ 22 ^] mnch weaker than 

the bonnds on K 2 k\) ref.|l], and we do not list them. Similarly, the existing 
bonnds on |A'{\^ ^' 22 k\^ l-^ 2 ifc ^'i 2 k\ mnch stronger than onrs. Fnrthermore with ttik- = 
(493.677 ± 0.016) MeV, rris = (122.5 ± 47.5) MeV and rris = (21 ± 4)md (see ref. [THj) we 
have the following new bonnds 
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( 20 ) 


The npper bonnd on Re[A)j\2 A212] obtained from two bonnds on the single conpling 
constants is stricter than the one we obtained. 


2.5 B- ^ if ^ 

For the charged R-meson decay the procednre is slightly different since it has not been 
directly measnred. Unlike the two previons cases one only has an experimental npper 
bonnd on the branching ratio R, see ref. jT^l, and thns has to go back to ea. (fTT|l . This 
has been done in ref. |I2|. We go beyond their work with a more conservative acconnt of 
the experimental errors and obtain weaker bonnds. We also work from the beginning in 
the mass eigenstate basis to avoid model dependent resnlts, see ref.jlH]. 

First / = 3. The theoretical predictions are limited by < 5.7 x 10“'^. 

As the total widths and T^ji^totai agree fairly well one has Ff^+S ^^i^totai^ so 
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that, utilizing eq. mj, we obtain for the branching ratio 
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To keep the combined uncertainties of iRisI and fs as small as possible we use the 
theoretical prediction (see ref. na) 
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In order to take into account the correlated uncertainties in V 13 /V 31 we use the Wolfenstein 
parameterization (see e.g. ref. my 


Ri3 _ p-iy 
Vsi 1 — ^ — p — iff 

The Wolfenstein parameters are given by (see ref. (201) P = 0.21 ±0.12, y 


(23) 
0.38 ± 


0.11, A = 0.222 ± 0.004, all at 95% C.L.. We thus obtain for the theoretical prediction 
/^^totai = (1-05 ± 0.65) X 10“^. The lower value should be used in eq. (j2T| . to 
be compared with the experimental upper bound. Thus 


2 Re[iA^ 3 i 33 ] ± E |-^313np < 13.3. 

n 


(24) 


In the following, we again assume that only two coupling constants are non-zero. Thus we 
have |i^ 3 i 3 i|, |R" 3132 | < a/ 13.3. Furthermore the imaginary part is again taken to be about 
the same as the real part, hence —± | < Re[ 7 ^ 3133 ] < \J\ + Thus, with 

IW 3 I = 0.0035±0.0015 (see ref.[ISl), ^b- = (5279.0±0.5) MeV, ruf, = (4200±200) MeV 
and nir = (1776.99 ± 0.29) MeV, see ref.[Tn] , we obtain 
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According to ref. [12] the bounds on |A'/i 3 Aifc 3 |, |A(j*ifc A; 3 /, A' 23 /, Re[A' 3 /fc A' 33 J and 

Re[A'i*i 3 A 313 ] are not better than the previous ones; furthermore the bound on |A '(\3 A 213 I 
is weaker than the product of the two bounds on the single coupling constants, see ref. | 3 |. 
Analogously, for / = 1 , 2 , 


(1 ± 2 Re[iF3i//] ± E l^3i/n/) X ^B- < 5.7 x 10“/ (26) 

n 

where tb- is the R-meson life time. Instead of arguing that the error on TE ,c is 
1 see the third line below ea. lO!^ . we are going to be as conservative as 
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possible. Due to isospin invariance and B~ have the same decay constant. From 
ref.j^I], /so = (200 ± 30) MeV, and thus with our convention /s = (141 ib 21) MeV {c.f. 
footnotelH). Therefore, ll/sp = (0.24 ± 0.22) MeV^, and with tb- = 1.655 x 10“^^ s 
(see ref.jin]), we obtain for / = 1 that (1 + ...)(9.0 ± 8.3) x 10“^^ < 1.5 x 10“^ and for 
/ = 2 that (1 + ...)(3.8 ± 3.5) x 10“^ < 2.1 x 10“®. Working with the lower value, for 
/ = 1 we get 
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The bounds on |A )./3 A 2 A:i| and Re[A )./3 Aifci] are not improved compared to the previous 
ones, see ref. [T2j, and the ones on |A'(\;i, A 23 fc|, |A'(\^ A 33 ;i,| and Re[A'(\^ A']^ 3 ;j] are too poor 
to be listed. / = 2 yields 
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Ref.jT2] states that there exist better bounds on |A'y'i 3 Aife 2 | and Re[A '/^3 A 2 A: 2 ]- The bounds 
on |A 2 ifc A)^ 3 ^|, |A 2 *ifc ^ 33 ^ 1 , Re[A 2 ifc A 23 fc] are almost the same as the single bound on |A 2 ifc|. 


3 ^p-Decay of Neutral Mesons 

3.1 Calculation of the Bounds 

Now we deal with the bound state decaying into and with momenta 

PitP 3 tP 2 : respectively. We consider only n ^ f, in which case the process does not 
occur in the SM and therefore no contributions from loop-diagrams have to be taken 
into account. We proceed as in the previous section. In eq. the 9*^ term together 
with its c.c. contributes to the decay, in analogy to the d/j^-exchange in the last section. 
Furthermore the 2”^ term together with the c.c. of the 5*^ and the 5*^ term together with 
the c.c. of the 2’^'^ contribute, both in analogy to the ^^^-exchange. The Hamiltonian is 
given by 
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Using the results corresponding to eq. o we obtain 
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Due to the large experimental error in we can neglect m^™ compared to m^/ 

(with f,n chosen correspondingly) and m^/ compared to Thus, focusing again 

on the bounds on products of two coupling constants, with all other coupling constants 
vanishing. 


I ^ijfn I) 


IH, 


ijfn I 


\B 


jinf\ 


mef 


mef 


< 


20 


f(dj^di) \ 


-pcxp. upper bound j j^exp 
{di^d*)^ef+endi / ^ {d7^d')-Atotal 


(34) 


Here r is the mean life time. The same considerations apply to mesons that have wave 
functions of the form 


<• = 


■B\jdi°dj ± (rfV) 


(35) 


one replaces every Aijfn by -^{Aijfn ± Ajifn), and likewise for Bijfn, B*inf- As in the 
previous section, we will apply eq. (ED only to processes with satisfactory experimental 
data, as was done in ref. dn, treating among other processes B^ ^ we conhrm 

their results. 
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3.2 B 
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We now consider > /i + e*", i, / = 2, j = 3, n = 1 and i j. The relevant 

parameters are given by/^o = (1.16±0.04) /^o, seeref.j2I], /i + e*^) < 6.1x10“®, 

Tfio = (1.464 ± 0.057) X 10“^^ s, m^o = (5369.6 ± 2.4) MeV, see ref.[T3]. Thus 
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(36) 


Our results for \ 2 k 2 \ 7^ 2) and fc = 1 for the second bound are weaker than the 

products of the bounds on single coupling constants, see ref.jl]. 
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77® is dehned as [Tf® + e77°]/Vl + with 77 ®2 = =*= 77®^]/\/2. e parameterizes the 

CP-violation. If we neglect e, 77® = [77® — 77®'"] /a/2, with 77® = {s'^d). From ref.[T3] one 
has = (497.672 ±0.031) MeV, = (5.17±0.04) x 10“® s and B{Kl ^ p±e®^) < 
4.7 X 10“^^. Ref. [13] gives fK = (159 ± 1.4 ± 0.44) MeV, which in the convention we use 
gives the central value 112.4 MeV. Hence, the hrst two bounds updating previous ones. 




/* I 

2fcll 


< 3 X 10 


-7 I 


m- 


Ul" 


100 GeV, 


7 * I 
'2fc2l 


< 3 X 10 


-7 


m- 


lA'i.i A 
1^21 Afci2| < 6 X 10“® 
|Al/2Afci2| < 6x10“® 
< 6 X 10“® 


Ul'" 


100 GeV, 


m- 


VL" 


100 GeV, 


m- 


vl" 


100 GeV, 


|Afci2 A 


m- 


VL 


A:2ll — 


11 


100 GeV 














W2I -^^211 — 


< 6x10“® I 


m- 


I'L 


100 GeV, 


(37) 


3.4 


TT 


/i 




With small modifications the result of section o can also be carried over to admixtures 
of with as the latter term does not contribute to any decay because the 

M-type quarks do not couple together to the operators. However, we shall limit 
ourselves to the 7r°: r] and ri' are more complicated, see ref.^H], and the experimental 
data do not suffice to extract satisfactory bounds. 

The relevant parameters here are 171^,0 = (134.9766 ±0.0006) MeV, 7^,0 = (8.4 ±0.6) x 
10“^'^ s and B{'k^ —>■ p ± e^) < 3.8 x 10“^°, see rebjlHI. Thus 
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,100 GeVy 

In ref.jlj a much stronger bound is stated for A 2 fci|. Furthermore they present a 
better bound for A^ 2 ilj from ref.jlj one hnds a stricter bound on IA 211 A 212 I, 
based on the bounds on single coupling constants. 


4 Summary 

We have determined the bounds on products of ftp coupling constants from leptonic 
meson decays. In many cases these bounds are better than previous bounds. We have 
summarized the bounds in the tables at the end of this text. With the formulae given the 
bounds can easily be updated when the data improve. Furthermore, if additional decays 
are measured (e.g. from the H-factories) one can determine additional bounds. Fq. |T7j) 
and eq. m can be used to consider 12 cases: D {i = 1, j = 2), {i = 2, j = 2), B 
{i = 3, j = 1), B~ {i = 3, j = 2) decaying into e ± and ii + {f = 1, g = 2), e + 
and T + {f = 1, g = 3), g. + and t + {f = 2, g = 3); eq. ljH^ can be applied to 
the decay of {i = 2, j = 3) to t + {f = 3,n = 1) or t + {f = 3,n = 2), and the 
decay of the T, (i = j = 3) to r ± e*" or r ± /i*" or p ± e'" (/ = 2, n = 1). 
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